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Theory predicts that the currents in scanning tunneling microscopy (STM) and the attractive 
forces measured in atomic force microscopy (AFM) are directly related. Atomic images obtained in 
an attractive AFM mode should therefore be redundant because they should be similar to STM. 
Here, we show that while the distance dependence of current and force is similar for graphite, 
const ant- height AFM- and STM images differ substantially depending on distance and bias voltage. 
We perform spectroscopy of the tunneling current, the frequency shift and the damping signal at 
high-symmetry lattice sites of the graphite (0001) surface. The dissipation signal is about twice as 
sensitive to distance as the frequency shift, explained by the Prandtl-Tomlinson model of atomic 
friction. 



The capability of scanning tunneling microscopy 
(STM) Q and atomic force microscopy (AFM) to 
resolve single atoms in real space makes them powerful 
tools for surface science and nanoscience. When operat- 
ing AFM in the repulsive mode, protrusions in the images 
simply relate to the atoms because of Pauli's exclusion 
law. In contrast, the interpretation of STM images is 
more complicated. The Tersoff-Hamann approximation 
0, valid for tips in an s-state, interprets STM images as 
a map of the charge density of the sample at the Fermi 
energy. Depending on the state of the tip, atoms can ei- 
ther be recorded as protrusions or holes, and tip changes 
can reverse the atomic contrast 0, . Theoretical pre- 
dictions regarding the relation of forces and tunneling 
currents / state that tunneling currents and attractive 
forces are directly related, thus AFM would not pro- 
vide any new physical insights over STM. Chen has 
found that the square of the attractive energy between tip 
and sample should be proportional to / with experimen- 
tal evidence in @|. Hofer and Fisher 7] suggested that 
the interaction energy and / should be directly propor- 
tional, experimentally found in H, 0. In this Letter, we 
investigate the experimental relationships between tun- 
neling currents and conservative as well as dissipative 
forces for graphite probed with a W tip by performing 
local spectroscopy on specific lattice sites. While force 
spectroscopy on specific lattice sites |l£| and combined 
force and tunneling spectroscopy on unspecific sample 
positions 0, Q have been performed before, the mea- 
surements reported here encompass site-specific spectra 
of currents and forces, supplemented by simultaneous 
constant-height measurements of currents and forces that 
allow a precise assessment of the validity of the theories 
regarding currents and forces in scanning probe experi- 
ments. 

In graphite (see Fig. 1), the electrons at Ep are 




FIG. 1: (color) Crystal structure of graphite. The unit cell 
(green) consists of two layers with inequivalent basis atoms a 
(white) and (3 (red). The a-atoms have direct neighbors in 
the adjacent atomic layers as indicated by the dotted lines, 
the /3-atoms are above a hollow site (h). (a) Perspective view, 
showing three layers formed by hexagonal rings, (b) Top view 
with surface unit vectors u and v. The line A = — v + u 
connects the a, (3 and h lattice points, spaced by 142 pm. 



concentrated at the (3 sites, and only these atoms are 
'seen' by STM at low-bias voltages. The state-of-the- 
art method for atomic resolution force microscopy is fre- 
quency modulation AFM (FMAFM) [13, 

where the fre- 
quency shift Af of an oscillating cantilever with stiffness 
eigenfrequency fo and oscillation amplitude A is used 
as the imaging signal 0, ^| . The bonding energy be- 
tween two adjacent graphite layers at distance a can be 
approximated by a Morse potential 
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with Ehond ~ —23 meV and k ~ 8nm _1 per atom 
pair 14]. The 'normalized frequency shift' ^(x^y^z) = 
kA 3 / 2 Af(x, y, z)j fo connects the physical observable Af 
and the underlying forces F ts with range A, where 7 ~ 
0.4F^A 5 (see Eqs. 35-41 in 13]). For covalent bonds, 
the typical bonding strength is on the order of -1 nN with 
A « 1 A, resulting in 7 w —4 fNyS, where a negative sign 
indicates attractive interaction. For graphite, the inter- 
layer bonds are much weaker and the potential of Eq. 1 
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FIG. 2: (color) Simultaneous records of a) tunneling current, 
b) frequency shift and c) damping. Image size 3 nm x 1 nm, 
tip bias 100 mV, scanning speed 40nm/s. 



results in 7 m i n = — 0.1fN^/m. The interaction of a tip 
atom with a graphite surface may be stronger than the in- 
terlayer bonds but should still result in — 7 m i n < 1 fN-y/m, 
posing a challenge for AFM imaging. True atomic reso- 
lution on graphite by AFM has so far only been obtained 
at low temperatures, first by Allers et al. using large- 
amplitude FM-AFM with Af = -63 Hz* f = 160 kHz, 
k 35N/m and A = 8.8 nm, thus 7 = -11.4 fNy/m. In 
spectroscopic measurements by the same group, a min- 
imum of 7 w —60 fN^/m has been observed (Fig. lb 
in [l£| ). Because this value is more than two orders of 
magnitude greater than the estimate above, it is expected 
that long-range forces have caused a large contribution in 
that experiment. Here, we use FMAFM with sub-nm am- 
plitudes which greatly reduces the influence of long-range 
forces and enables simultaneous STM operation [13 . 

We use a low-temperature STM/ AFM operating at 
4.9 K in ultrahigh vacuum 0. The microscope uses 
a qPlus sensor [l£j for simultaneous STM/ AFM oper- 
ation (k = 1800 N/m, f = 11851.75 Hz, quality factor 
Q = 20000). All data (spectroscopy and images) are 
recorded at A = 0.25 nm. The tip is prepared by dc 
etching (3 V) of a polycrystalline tungsten wire. The 
frequency shift is measured with a commercial phase- 
locked- loop detector (EasyPLL by Nanosurf AG, Liestal, 
Switzerland). The instrument is thermally well con- 
nected to a liquid He bath cryostat at 4.2 K, leading to a 
drift rate of « 20pm/hour. Non-conservative tip-sample 
interactions lead to damping, and the energy AE ts that 
has to be provided for each oscillation cycle to keep A 
constant is recorded simultaneously with / and Af. 

To check if the three data channels /, Af and AE ts 
are produced by the same tip atom, we have scanned an 
area that contains a step edge (Fig. |2j). The step edge 
appears at the same position in all three data channels, 
thus the signals are produced by the same tip atom. 

The physics of the interaction is best explored by 
performing /, 7 and A£^ s -spectroscopy at the high- 
symmetry lattice sites. Figure 01 shows I(z), 7(2) and 
AE ts (z) taken at the a-, /3- and h-lattice sites. All three 
signals initially increase roughly exponentially, as shown 
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FIG. 3: (color) Experimental spectra of /, 7 and AEts as a 
function of distance for a tip bias of 160 mV. The insets are 
logarithmic plots of /, —7 and AE ts - The grey curves are 
taken at position '1' in Fig. 4a), the black curves at position 
'2' and the red curves at position '3'. The insets are views of 
I(z), 7(2:) and AEts(z) for — 0.1 nm< z < on logarithmic 
scales. 



in the log-scale insets. The correspondence between ex- 
perimental images and the lattice sites is established by 
the analysis below. The current increases exponentially 
in the distance regime from z = to z = —100 pm, fol- 
lowed by a step-like increase for smaller distances. We 
assume that the tip is essentially in contact to the upper 
graphite layer for distances smaller than -100 pm. Given 
that the electrical conductivity of graphite is small at low 
temperatures and small in ^-direction, we conclude that 
the step-like increase in current for distances smaller than 
-100 pm is caused by an increasing number of graphite 
layers becoming available for charge transport. The cur- 
rent spectra in Fig. 3a) are maximal for most distances 
on site '3'. Because site '3' is a current maximum, we 
identify it as a /3 position. 

The normalized frequency 7 shift decreases down to 
a distance of ~ —170 pm, followed by a slight increase 
down to « —550 pm and a sharper increase for even 
smaller distances. For z < —800 pm, 7 remains con- 
stant because the cantilever remains in contact to the 
sample during the entire oscillation cycle. In the z-range 
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from to -100 pm where the short-range chemical bond- 
ing forces start to emerge, the magnitude of 7 increases 
exponentially (see inset). It is believed, that the short- 
range forces between AFM tips and graphite originate 
from van-der-Waals forces [l^ with their typical 1/z 7 - 
distance dependence. The experimental data shows that 
when the interatomic distance approaches the atomic di- 
ameters, an exponential force dependence prevails. 

The damping signal initially also increases exponen- 
tially for z < 0, reaches a plateau for z < —170 pm, 
decays to zero from z < —600 pm and remains zero for 
z < —800 pm because the cantilever remains in contact 
for the whole oscillation cycle. This points to a damp- 
ing mechanism as described by Prandtl [2fl[ and Tom- 
linson |2l|, where the energy loss is caused by a pluck- 
ing action of the atoms on each other. The energy loss 
per cycle is simply related to the maximal attractive 
force F ts min and the stiffness of the sample k sarnp i e with 
AE ts =Fl mi J(2k sample ) 0. 

The insets are logarithmic plots of I(z), 7(2) and 
AE ts (z) for — 0.1 nm< z < 0, showing an almost expo- 
nential distance dependence in that range. For the tun- 
neling current, we find a decay constant kj = 13 nm -1 
at the /3-site and kj = 15 nm -1 at sites 'l'and '2', lead- 
ing to an apparent barrier height of ~ 2eV. The decay 
constant of 7 and thus the interaction potential (Eq. 39 
in 0) is ft 7 = 12 nm -1 at the /3-site and ft 7 = 16 nm -1 
at sites 'l'and '2'. The decay constants for kj and ft 7 are 
equal within the measurement accuracy, thus the theory 
by Hofer and Fisher Q appears to hold for the interac- 
tion of W with graphite. The damping signal decays with 
KAE ts — 20 nm -1 at the /3-site and kae ±s — 30 nm -1 at 
sites 'l'and '2' as expected from an energy loss propor- 
tional to the square of the attractive force. 

The normalized frequency spectra shown in Fig. 3b) 
are rather similar, expect for the grey curve recorded 
at site '1'. In the distance regime from —600pm to 
—800 pm, the grey curve is shifted by « —40 pm. If the C 
sample atoms and the W tip atom are assumed to be hard 
spheres with a diameter of 142 pm and 273 pm respec- 
tively, the W tip atom could protrude 56 pm deeper on 
top of the hollow sites than at a or (3 sites. We therefore 
conclude that position '1' (grey curves) corresponds to a 
hollow site (h in Fig. 1(b)). Interestingly, Fig. 3c) shows 
that the dissipation is significantly larger on the hollow 
site than on top of a or (3 sites. The three data chan- 
nels were acquired simultaneously and the range from 
z = — 1.7nm to 0.6 nm and back to — 1.7nm was ramped 
within 60 s. The sequence of the three sites was scanned 
three times, so a total of 6 spectra was collected for each 
a-, (3- and h-site. While the drift rate of our instrument 
is very low, piezo creep caused z-offsets of consecutive 
scans. These offsets were calibrated by comparison with 
constant-height scans which provide precise cross refer- 
ences for /, 7 and AE ts at the a-, (3- and h-sites for a 
given z-value. The local spectra at high-symmetry sites 
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FIG. 4: (color) Constant- height measurements of /, 7 and 
AE ts in attractive and repulsive distance regimes for a tip 
bias of 160 mV (a-c) and -60 mV (d-f). The hexagons show 
the proposed positions of a (white) and f3 (red) atoms. The 
brightness is proportional to |/|, 7 and AE ts . The inset in c) 
left shows a higher harmonic image [2^, indicating that the 
tip state is not perfectly symmetric with respect to the z-axis. 



were supplemented by constant height scans at various 
^-positions shown in Fig. 4. Figures 4a)-c) show /, 7 
and AE ts at a tip bias of 160 mV in the fully attrac- 
tive mode at z « —100pm (a), in a weakly repulsive 
regime at z « —350 pm (b) and a fully repulsive mode 
at z « —750 pm (c). If the attractive interaction be- 
tween tip and sample was only mediated by the electronic 
states that contribute to the tunneling current, the I- 
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and 7-images in Fig. 4a) should be exactly inverse. Evi- 
dently, this is not the case. We therefore conclude, that 
electronic states that do not contribute to the tunnel- 
ing current may contribute to attractive interaction. In 
the repulsive regime shown in Fig. 4c), the repulsion 
is strongest above the a sites, almost as strong on the 
(3 sites, and weakest on hollow sites, in agreement with 
with Fig. 3a) and the Pauli exclusion law. The current 
has a local maximum on top of the (3 sites, and AE ts 
has a pronounced maximum at the hollow sites. Because 
of the long time scales involved in damping measure- 
ments, atomically resolved energy loss measurements are 
prone to lateral shifts [24( for fast scanning. The scan- 
ning speeds used in the data of Fig. 4 were « 0.3nm/s, 
therefore time delays in the acquisition channels are neg- 
ligible and the /, 7 and AE ts images match precisely in 
forward and backward scans. In our previous simultane- 
ous STM/AFM measurement on graphite [Tjj, we found 
a lateral shift of 68 pm of the current maxima with re- 
spect to the corresponding 7 maxima. While the tip is 
not perfectly symmetric with respect to the z-axis (see 
caption Fig. 4), it is evidently more symmetric than in 
[13. Arai and Tomitori 25] have recently shown that 
force interactions are also a function of bias. Figure 4d)- 
f) shows constant-height images at a different tip bias of 
-60 mV. In the attractive regime shown in Fig. 4d), the 
/- and 7-images are approximately inverse, i.e. the local 
minimum in 7 coincides with the local maxima in /. In 
the repulsive regime shown in Fig. 4f), again the repul- 
sion is strongest above the a sites, almost as strong on 
the f3 sites and weakest on hollow sites. 

While the distance dependencies of current, force and 
dissipation as revealed by constant-height images and lo- 
cal spectroscopic measurements are qualitatively similar, 
the contrast observed in the const ant- height images is 
larger than expected from the local spectra. The reason 
for these subtle discrepancies is revealed by the constant 
height data shown in Fig. 4. Frequency shift and current 
images are different in all distance regimes, and slight 
shifts between the high- symmetry points in / and 7 im- 
ages are present. 

In summary, the spectroscopy experiments show that 
the ^-dependence of force and current is roughly the same 
as predicted in Hofer and Fisher's theory [3 for graphite. 
However, the constant-height experiments proof that at- 
tractive forces and currents are not directly related and 
STM and AFM do provide different information. The 
dissipation measurements reveal that the theories on 
atomic friction introduced by Prandtl [2£| and Tomlinson 
[2~H are the key mechanism for damping on the atomic 
scale when imaging soft samples. The spatial resolution 
that is possible by scanning probe microscopy scales with 
the decay lengths of a physical observables [2fj. At the 
onset of damping, the decay length of dissipation is only 
half the value of the force. This offers an alternate expla- 
nation of the impressive resolution obtained in damping 
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